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Abstract 

Cohomology operations (including the cohomology ring) of a geometric 
object are finer algebraic invariants than the homology of it. In the litera- 
ture, there exist various algorithms for computing the homology groups of 
simplicial complexes ([Mun84], [DE95, ELZOO], [DG98]), but concerning 
the algorithmic treatment of cohomology operations, very little is known. 
In this paper, we establish a version of the incremental algorithm for com- 
puting homology given in [ELZOO], which saves algebraic information, al- 
lowing us the computation of the cup product and the effective evaluation 
of the primary and secondary cohomology operations on the cohomology 
of a finite simplicial complex. The efficient combinatorial descriptions at 
cochain level of cohomology operations developed in [GR99, GR99a] are 
essential ingredients in our method. We study the computational com- 
plexity of these processes and a program in Mathematica for cohomology 
computations is presented. 

1 Introduction 

A simplicial complex is a well-known discrete model of a geometric object, which 
consists of a collection of simplices that fit together in a natural way to form 
the object. In order to classify simplicial complexes from a topological point of 
view, a first algebraic invariant that can be used is homology, which in some 
sense, counts the number of holes of the object. 

We can cite two relevant algorithms for computing homology groups H*K of 
a simplicial complex K in R": (1) the classical algorithm based on reducing cer- 
tain matrices to their Smith normal form [Mun84]; (2) the incremental algorithm 
[DE95, ELZOO, EZ01], avoiding the severe computational costs of the reduction 
to Smith normal form and consisting of assembling the complex simplex by 
simplex and at each step updates the Betti numbers of the current complex. 
Starting with the boundary of a negative simplex, this persistence process finds 
the cycle which is destroyed by this simplex through the search, computing in 
this way the geometric realization of a homology cycle. It runs in time at most 
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0(to 3 ), where m is the number of simpliccs of the complex. For simplicial com- 
plexes embedding in R 3 , this complexity is reduced to 0{m) in time and space 
[DE95]. The algorithm proposed in [DG98] is based on simulating a thicken- 
ing of a given complex in R 3 to a topological 3-manifold homotopic to it, and 
computing the homology groups of the last one using classical results. The time 
and space complexity is linear and this method also produces representations of 
generators of the homology groups. 

In general, computing homology is not enough for determining whether two 
geometric objects are homeomorphic or not. Finer algebraic invariants such as 
the cohomology (an algebraic dual notion to homology), the cup product on 
cohomology or cohomology operations [Spa81], allow us to topologically distin- 
guish two geometric objects having isomorphic homology groups. For example, 
a torus and the wedge product of a sphere and two circles have the same homol- 
ogy but the respective cup products on cohomology are "essentially" different. 
Using a field as the coefficient group, for example, Z2, the cohomology H*K of a 
simplicial complex K gives us the same topological information as the homology 
of it. However, the additional ring structure on the cohomology determined by 
the cup product and cohomology operations cannot directly be produced from 
the algorithms previously mentioned for computing the homology. Roughly 
speaking, a cohomology operation 9 : H m (—:G) — > H n (—;G') is a homomor- 
phism that acts on cohomology (G and G' being groups); relevant examples 
of cohomology operations are Steenrod squares, Steenrod reduced powers and 
Adem secondary cohomology operations [MT68] . As an example of the strong 
constraints that these operations impose on the cohomology of spaces, we can 
cite that the use of this machinery is essential for showing that there do not 
exist spaces X having cohomology H*(X;Z) a polynomial ring Z[a) unless a 
has dimension 2 or 4. 

In this paper, we make use of an explicit chain contraction (a special chain 
equivalence) connecting the chain complex C*K, canonically associated to a 
simplicial complex K and its homology H^K. Moreover, from this datum we can 
derive a cochain contraction from the cochain complex C* K = Hom(C^K; Z2), 
to the cohomology H*K. Using this information, we can compute: 

1. Geometric realizations of (co)homology generators. 

2. The (co)homology class of a (co)cycle in terms of (co)homology generators. 

3. The construction of a (co)boundary of a given (co)cycle. 

4. The induced homomorphism at (co)homology level of a simplicial map 
between two complexes. 

5. The cup product on cohomology and some primary and secondary coho- 
mology operations. 
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The first problem is to construct such chain contractions from C*K to H^K. 
In [GR01], a translation of the classical matrix algorithm (1) in terms of chain 
contractions is designed. In this paper, we design a version of the incremental 
method described in [ELZOO] in terms of chain contractions. The complexity 
of our method is also 0(m 3 ) where m is the number of simplices of K, but 
our algorithm saves information which allows us, for example, to compute the 
following operations: 

1. The cohomology ring of K in 0(m ). 

2. The Steenrod square operation Sq l a n of a cohomology class a n of degree 
n in 0(i n - l+1 m) (see [GR99a]) 

3. The Adem secondary cohomology operation ^> 2 a2 of a cohomology class 
a 2 G KerSq 2 H 1 (K;Z 2 ) in 0(m 3 ). 

In fact, the modus operandi for evaluating a mod 2 cohomology operation O : 
H m K — > H n K on a cohomology class a m is the following: 

1. First, given a finite simplicial complex K, construct the chain contrac- 
tion from C*K to H*K (denoted {f*,g*,(j)*) : C*K => H*K), using our 
version of the incremental technique. 

2. Evaluate O on the cohomology class a m using the diagram 

C*K £ H*K 

C*K 4 H*K, 

where O : C* K — >• C*K is a cochain operation associated to O whose for- 
mulation is explicitly given in simplicial terms. An efficient combinatorial 
description O for O being a Steenrod square [GR99, GR99a], a Steenrod 
reduced power [GR99] or some Adem secondary cohomology operations 
[GR01] have already been done by the authors. We do not deal with 
this question in this paper, but it is necessary to say that the algorith- 
mic approach we give here will only be valid if combinatorial pictures of 
cohomology operations at cochain level are determined. 

Let us observe that in this paper we deal with the general case of R™. 
Versions in terms of chain contractions of the algorithms given in [DE95] and 
[DG98], designed for the special case of R 3 , would allow us to considerably 
reduce the computational costs of the processes. 

2 Homology and Chain Contractions 

In this section, we design a version of the incremental algorithm of [ELZOO] 
in terms of chain contractions. In this way, we construct a chain contraction 
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from the chain complex canonically associated to a simplicial complex K, to its 
homology. Let us observe that passing to cohomology is not a problem if we 
use a field as the ground ring. The resulting cochain contraction from C*K to 
H*K will help us to compute the cup product on cohomology and cohomology 
operations. 

Now, we give a brief summary of concepts and notations. The terminology 
follows Munkres [Mun84] . 

Throughout this paper, we consider Z 2 is the ground ring and \x denotes 
the product on Z 2 . A q— simplex a in R" (where q < n) is the convex hull of 
q+l affincly independent points {vo, v q }. We denote a = (vo,...,v q ). The 
dimension of a is \a\ — q. A 0-simplex is a vertex, a 1-simplex is an edge, a 
2-simplcx is a triangle, a 3-simplex is a tetrahedron, and so on. An i-face of 
a = (vq, v q ) (i < q) is an i-simplex whose vertices are in the set {«o, • v q }. 
The (q — l)-faces of a are called the facets of a. A simplex is shared if it is a 
face of more than one simplex. Otherwise, the simplex is free if it belongs to 
one higher-dimensional simplex, and maximal if it does not belong to any. A 
simplicial complex if is a collection of simplices such that: 

• If r is a face of a G K, then t E K. 

• If a', a E K, then a' D a E K or a' n a = 0. 

Let us notice that K can be given by the set of its maximal simplices. The 
dimension of K is dimK = max{\o-\ : a G K}. In this paper, all the simplices 
have finite dimension and all the simplicial complexes are finite collections. The 
set of all the g-simplices of K is denoted by K^ q \ If L is a subcollection of 
K that contains all faces of its elements, then L is a simplicial complex in its 
own right; it is called a subcomplex of K. Let K and K' be two simplicial 
complexes. A map / : — > K'^ such that whenever (w , —iVq) G K then 
f(vo), f(v n ) are vertices of a simplex of K' , is called a vertex map. 

Algebraic Topology is the study of algebraic objects attached to topological 
spaces; the algebraic invariants reflect some of the topological structure of the 
spaces. 

The chain complex C*K associated to a simplicial complex if is a family 
{C q K, d q } q >o defined in each dimension q by: 

• C q K is the free abelian group generated by the q-simplices of K. An 
element a = o\ + ■ ■ ■ + a m of C q K ioi G K^) is called a q-chain. 

• d q : C q K — > C q -\K called the boundary operator is given by 



where (vo,—,v g ) is a g-simplex of K and the hat means that Vi is omitted. 
By linearity, d q can be extended to C q K, where it is a homomorphism. 
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A g-chain a is called a q-cycle if da = 0. If a = db for some b £ C q +\K then 
a is called a q-boundary. We denote the groups of g-cycles and g-boundaries 
by Z q K and B q K respectively, and define ZqK = CqK . Since B q K C Z g .ff, we 
can define the gi/i homology group to be the quotient group Z q K/ B q K , denoted 
by H q K. Given that elements of this group are cosets of the form a + B q K, 
where a £ Z q K, we say that the coset a + B q K, denoted by [a], is the homology 
class in H q K determined by a or a is a representative cycle of [a]. Let K and 
L be two simplicial complexes. A chain map f : C*K — > C*L is a family of 
homomorphisms 

{fq ■ C q K ->■ Cqi}q>0 

such that c^/q = f q -\d q for all g. Observe that for every vertex map f : K^ — > 
, we can obtain the corresponding chain map /# : C*K — > C*L such that 



U(vo,..,v q ) = { W*>'~' 



{f{vo),-,f(v q )) if f(vi) distinct 
otherwise 



Let h and k be two chain maps from C*K to C*L. A ckin homotopy from 
/i to fc is a family of homomorphisms 

{(f> q : C q K -> C 9+ iL} g >o 

such that d q +i<p q + 4> q -id q = h q + k q . Wc write ft, ~ fe if a chain homotopy 
between h and fc exists. Two chain complexes C*K and C*L are chain equivalent 
if there exist two chain maps / : C*K — > C*L and g : C*L C-*K such that 



/ff ~ lo.i and 5/ - 1 



C. K ■ 



Observe that, in this case, 4> q : C q K — > C q +\K for all q > 0. A ckm contraction 
[EM52] from C*K to C*Z/ is a chain equivalence such that 

fg = 1o»l and g/ - 1 C , K (that is, 1 C „ K + 5/ = 90 + 09) 

and 4> has the following "annihilation" properties: f(f> = 0, (f>g = and 00 = . 
We denote such chain contraction as (f,g,(f>) ■ C*K => C*L. Observe that if a 
chain contraction from C*K to C*L exists then L has fewer or the same number 
of simplices than K. Now, we show some examples of contractions. 

(a) Edge Contractions. 

Conditions under which edge contractions are homeomorphisms appear in 
[DEGN99]. Here, we show one condition under which edge contractions 
become, at algebraic level, chain contractions. 

Let if be a simplicial complex and r = (a, b) an edge in K. An edge 
contraction is given by the vertex map / : K^ — > = K^ — {a, b} U 
{c} where /(a) = /(&) = c, and f(v) — v for all v ^ a,b. 
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Let B be a subset of K that is not necessarily a subcomplex. Define 



B = {a'EK: a'<aEB}, St B = {a G K : a > a' G B} , 



Lk B — St B - St B , 
where a' < a means that cr' is a face of cr. 

If Lka n Lfc6 = Lkr, then a chain contraction (f#,g,4>) from C*i4T to 
C*L is defined as follows: 

— /# is the chain map induced by the vertex map /. 

— g : C*L — > C*K is such that 

gr = t \/t £ Stc, 
5(c) = (a), 

ujU (a) if we Lka, 

ujU (b) +wU (a, 6) if ueLkb-Lkr 



g(uj U (c» 



ui E Lkr and a) < w, 
w U (6) if uj E Lkb - Lk t, 

flu) < lu and ui E Lkr. 



- <f> : C^K -> C* + iif is given by 

0(u o , ...,v q ,b) = (v , ...,v q ,a,b) if (u , • v 9 ) £ Lkr 
and </>t = otherwise. 

(b) Simplicial Collapses. 

Suppose K is a simplicial complex, a E K is a maximal g-simplex and a' 
is a free (g — l)-face of a. Then, K simplicially collapses onto K — {a, a'}. 
More generally, a simplicial collapse is any sequence of such operations. A 
thinned simplicial complex M sco \{K) is a subcomplex of K with the con- 
dition that all the faces of the maximal simplices of M sco \(K) are shared. 
Then, it is obvious that it is no longer possible to collapse. There is an 
explicit chain contraction from C*K onto C* (M sco \K) [For99]. The follow- 
ing algorithm computes M sco \K and the chain contraction from C*K onto 
C*(M sco iK). Suppose that K is given by the set of its maximal simplices. 



Initially, M sco \K = K , tfrr = , fr = gr = T for each teK. 
While there exists a maximal simplex a with a free face a' do 

M scol K = M scol K - {a, cr'}, 
4>a' = a , fa' = cr' + dcr and fa = 
End while 
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(c) Contraction to a Vertex. 

Let a — (vo, ■ ■ ■ , v q ) be a simplex and let K[a] be the simplicial complex 
whose maximal simplex is a. It is obvious that we can obtain a chain 
contraction from C*K[<j] to (v ) using simplicial collapses. But now, we 
show another contraction from C*K[a] to {v ) determining the acyclicity 
of the simplex a. This last chain contraction is the key for constructing 
another one from any simplicial complex to its homology as we will see in 
the following section. We dehne (f a ,g<j,4>a) '■ C*K[a] => (v ) as follows: 

f<j( v i) — i v o) 0<i<g, and / ct (t) = otherwise; 
9a(vo) = (vo); 

<j} a {vQ,v h , . . . ,v jn ) = and 4> a {v h , . . . , v jn ) = (v , v j± , . . . , v jn ) 
where 1 < ji < • • • < j n < q. 

Let us observe that in this case (v ) represents the unique class of homol- 
ogy in H*K[<j]. 



2.1 Incremental Homology Algorithm and Chain Contrac- 
tions 

Our algorithm for computing a chain contraction from the chain complex of a 
simplicial complex K to its homology is based on the incremental algorithm for 
computing the persistence of the Betti numbers developed in [ELZOO]. 

The input of our algorithm implemented in Mathematica is the sorted set 
of all the simplices, K = {cti, . . . , <r m }, with the property that any subset 
of it, {ci, . . . , CTj}, i < to, is a simplicial complex itself. The output t = 
contraction IK] is a list of sorted lists. Each sorted list has three elements. 
The first one is a simplex a of K , the second one is the image of a under / 
and the third one consists of the image of a under <j). We omit in the list the 
simplices such that the image of them are null under / and <j>. In general, a 
class of homology a is represented by a simplex r, so in order to obtain the 
image of a under g, we only have to compute a — r + (j>dr. Moreover, a will be 
a representative cycle of a. 

Now, let us suppose we have constructed the list I ^contraction [L] for 
L = {<7i, . . . ,<7j_i}, i < m (if L = 0, we assume I = 0). We construct 
contraction [{<Ti, . . . , Uj}] as follows: 

If f [8^,11=0 then, 

£U {(<7j,<7j,(/>0-j)}, 

Else 

Replace [ Replace [ £, 

Solve [fW<Ti,t]=0] 

], 
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Solve [0 \_d(Ji, £~] =£7,] 

] 

End if 

where, for a simplex r, f[T,£] and 0[t, are, respectively, the second and 
the third clement of the list of £ that has r as the first element. If this list 
does not exist, then /[t, £1— and </>[t, £]= 0. Now, let us explain what 
contraction [{(Ti, Ci}] computes. If f[dai,£~\= then oi "creates a cy- 
cle", so in fact, m is a new generator of homology. Otherwise, f[d<Ji,£~\ is a 
sum of elements of the form J2a -eNcL °r The idea °f ^his last case is that Oi 
destroys the cycle generated by d<Ji in L. Therefore, we impose /[9tTi,£]= 
and (j) [d<Ji, £]= <Ji. We replace these relations in £ with the commands Replace 
and Solve. 

At the end of the algorithm, all the elements of the form <pr are replaced 
by zero. For obtaining the morphism g and the representative cycles of the 
homology classes of K, we compute r + (j)dr for each simplex r (the generators 
of homology) satisfying that /[r, £] = r in the list £ = contract ion [K\. We 
create a new list of sorted lists, called representativeCycles[if ] such that in 
each sorted list the first element is a generator of homology, r, and the second 
element is its image under g, r + </><9t. Observe that this last chain is, in fact, 
a cycle: 

<9(t + (f>dr) = dr + d(f>dr 

= dr + (gf - 1 - 4>d)dT 

= gfdr — <pddr [ since ddr — 0, then ] 

= gfdr [ since, by construction, fdr = 0, then ] 

= 0. 

It is easy to check that (f 7 g 7 (f>) is, in fact, a chain contraction from C*K to 
H-xK. Observe that given a cycle a, if fa = then a is also a boundary. In 
order to compute a chain a' such that a = da', we can use the relation 

a ~ 9f a = 4>da + d(f>a . 
Since da — and fa = 0, we have a = d<j)a . 

Theorem 3 The complexity of our algorithm for computing the homology of 
a finite simplicial complex K and a chain contraction from C^K on H*K is 
0(m 3 ), where m is the number of simplices of K. 

Proof. 

Let K = {<7i, . . . , a m } and d = dimK. Suppose that we have computed 
£ =contraction[{<7i, . . . , . In the worst case, we have to solve / [cVj, £1 

= and (f> ld<Ji , £~\ = &i. Observe that the number of simplices involved in d<7i is 
less or equal than the dimension of Oi which is at most d and then, the number of 
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simplices involved in the formulas of / \.doi,l~\ and <j) Ldai, £1 is 0(dm) = 0(m). 
Since we have to solve the equations and replace the solution in £, the total cost 
of these operations is 0(m 2 ). Moreover, for obtaining the representative cycles, 
we have to compute r + 4>dr for every generator of homology. The cost of this 
is also 0(m 2 ). Therefore, the total algorithm runs in time at most 0(m 3 ). □ 

4 Cohomology and Cohomology Operations 

One reason in order to use the cohomology for distinguishing spaces instead of 
homology, is that the cohomology has additional structures, such as cup product 
and cohomology operations. If two spaces have isomorphic (co)homology groups 
but the behaviour of the ring structure or cohomology operations is different, 
then they are not homeomorphic. In this section we explain how we can compute 
the cup product and cohomology operations starting from a chain contraction 
from an algebraic object to its homology. We first need to define more concepts. 
The cochain complex associated to K, denoted by C*K, is the family 

{CK,5*} q > , 

defined in each dimension q by: 

• The group C q K = Hom(C q K; Ti2)—{c : C q K — > Z2 is a homomorphism}. 

• The homomorphism S q : C q K — > C q+1 K called the coboundary operator 
given by 

S q ca = cd q+ \a 
where c G C q K and a 6 C q+ \K. 

The elements of C q K are called q—cochains. Observe that a g-cochain can be 
defined on and it is naturally extended by linearity on C q K. Z q K and 
B q K are the kernel of S q and the image of <5 9_1 , respectively. The elements in 
Z q K are called q-cocycles and those in B q K are called q-coboundaries. The qth 
cohomology group 

H q K = Z q K/B q K 

can be defined for each integer q. Take into account that since the ground ring 
is a field, the homology and cohomology of K are isomorphic. Moreover, given 
a generator of homology, a, of dimension q, we can define the corresponding 
generator of cohomology a* : H q K — > Z 2 such as 

a*a = 1 and a* j3 = for a ^ ft E H q K . 

One can also define the dual concept of chain maps and chain contractions, in 
the obvious way. Furthermore, starting from a chain contraction (f,g,(i>) from 
C*K to H*K, we construct a cochain contraction {f*,g*,4>*) from C*K to H*K 
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as follows. Let c G C*K and a* G H*K. Define f*c = eg, g*a* = a* f and 
4>*c = ccj). 

The cohomology of K is a ring with the cup product 

— : H l K®WK ^ H l+1 K 

defined at a cocycle level by (c ' c')<r = n(c(v , . . . , Vi)®c' (vi, . . . , where 
c and c' are an i-cocycle and a j -cocycle, respectively, and a = (vq, . . . , Uj+j) G 
(*+.?') i s suc h that «o < • • • < v i+ j. Using the chain contraction (f,g,4>) from 
C*K to H*K, we can compute the cohomology ring of K in the following way: 

Take a* and j3* , cohomology classes of K 
For every 7 G Hi + jK 

compute ((a*/) - (P* f))n 
End for 

Notice that the resulting cohomology class is determined by the cocycle c = 
(a*f) w (/?*/). 

In order to compute a cohomology operation O : H*K — > H* +t K, on one 
hand, we need to compute contraction [if] in order to obtain a chain con- 
traction (f,g,4>) from C*K to its homology and, on the other hand, we need 
a simplicial version O : C*K — > C* +l K of O. Therefore, for obtaining <D(a*), 
where a* G H*K 1 we only need to compute 0(a*f)g (for more details, see 
[GR01]). For example, from the combinatorial formulae of Steenrod squares 
given in [Ste47, SE62], 

Sq i : H*K ->• H* +l K , 

for calculating the cohomology class Sq l {a*) with a* in H q K, we only have to 
compute Sq l (a* f)g. More concretely, at cochain level, Sq t c = c ^ q -i c mod 2. 
Moreover, given a p-cochain c and a g-cochain c', c — „ d is a (p+q — n)-cochain 
defined by 

(c—„c')<7= M C ( U j cvcnZ 3 ) ®c'(Uj odd^)) 

0<io<"-<*Ti<p+Q — n 

where ct = (u , • • • , Up+ g -»), «o < ■•■ < Wp+g-i! ^° = (u , •■• ,«»<,), ^ = 
• • • ,%), for 1 < j < n, and = (v in , . . . ,v p+q - n ). Finally, we 

can express Steerond squares in a matrix form due to the fact that these coho- 
mology operations are homomorphisms. The process of diagonalization of such 
matrices can give us detailed information about the kernel and image of these 
cohomology operations. This information will be very useful in the next section 
in order to compute Adem secondary cohomology operations. 

5 Adem Secondary Cohomology Operations 

For attacking the computation of secondary cohomology operations, we will see 
in this section that the homotopy operator </> of the chain contraction (f,g,^>) 
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from C*K to the homology of K, is essential. 

First of all, we will need the following mod 2 relation [Ste47]: 

<5(c -„ c) = c -~-„_i c + c -~-„_i c + <5c — -„ c' + c w n (5c (1) 

where c and c' are two cochains. Now, we shall indicate how Adem secondary 
cohomology operations 

: N q K -> H q+3 (K; Z 2 )/Sq 2 H q+1 (K; Z), g > 2 

can be constructed (see [Ade52, Adc58]). N q K denotes the kernel of Sq 2 : 
H q (K; Z) — > H q+2 (K; Z2). These operations appear using the known relation: 

Sq 2 Sq 2 a + Sq 3 Sq 1 a = 

for any a <E H*(K; Z). For this particular relation there exist cochain mappings 

Ej : C*(K x K x K x K) ->• C*"^ 

such that mod 2 

(C — 9 _ 2 C) (C — g-2 C) + (C — g _l c) (c c) = (5£ 3(Z _3C 4 , 

where c is a g-cochain with integer coefficients. Making use of the relation (1) 
we have that mod 2 

(c — g -2 c) — 9 (c — 9 _ 2 c) = <5(fo — q 5b + b — q _! b) 

(C — g_i C) — ,_2 (C C) = 6(r) Sfj + 1] V) 

where b is a (q + l)-cochain such that c ^ q -i c = Sb and 77 = |(c ^ 9 c + c). 
Therefore 

{^3<J-3C 4 + & b + b —'g Sb + T) ^g-2 8r) + T) ~~^g-3 V, 1 > 2 

E 3 c 4 + b^ 1 b + b^ 2 5b + r)^5ri, q = 2 

is a mod 2 cocycle. If c is a representative g-cocycle of a cohomology class 
a G A^if with integer coefficients then, 

«I> q a= [w] + Sq 2 H q+1 K . 

Now, suppose Z 2 is the ground ring and suppose we have computed the con- 
traction (/, g, <j>) from C^K to H*K, £ ^contraction [K~\ . Then, the cochain b 
is <fi*(c -—-q-2 c) = (c w 9 -2 c)cp. Observe that for computing * 9 a*, a* G H*K 1 
we need to have a combinatorial expression of the morphism E 3q _ 3 . A method 
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for obtaining "economical" combinatorial formulae for E 3q - 3 is given in [GonOO]. 
For example, 

(£; 3 C 4 )(T = fJ,(c{Vo,V2,V 3 )®c(Vo,V 1 ,V2)®c{V3,V4,V 5 )®c{V2,V3,V 5 ) 

+c(v ,v 4 ,v 5 ) ® c(v 3 ,v 4 ,v 5 ) ® c(v ,v 1 ,v 2 ) ® c(v ,v 1 ,v 2 ) 
+c(v ,v 1 ,v 5 ) ®c(v 3 ,v 4 ,v 5 ) ®c{vi,V2,v 3 ) ®c(vi,V2,v 3 ) 

+ c(v 0l Vi,V 2 } ® c(v 2 , Vi, V 5 ) ® c(v 2 ,V3,Vi) ® c{v 2 ,v 3 ,Vi) 

+c(v 0l vi,v 2 ) ® c(v 2 ,v 3 ,v 5 ) ® c{v3,Vi,v$) ® c{v3,va,v$)) , 

where c is a 2-cochain and a = (vo, i>i, v 2 , V3, V4, v$) is a 5-simplex such that 
v < v\ < v 2 < v 3 < V4 < V5. Therefore, the steps for computing ty q are the 
following: 

1. Take a* G N q K making use of the diagonalization of the matrix of 
Sq 2 H q K. 

2. Compute c — a*f. 

3. Compute b = (c c)<f>, r\ = h(c ^ q c + c), b ^ q -i b, b -~^ q 5b, 
V V, V St/ and E 3q _ 3 c 4 . 

2. Compute wg. 

Let us explain with more detail the first step. In our implementation in Math- 
ematica, the command hclass[£, q] computes the list of all the cohomology 
classes of K in dimension q. We compute Sq 2 a* for each a* G hclass [£, q~\ and 
we write the result as a vector sq2 [ £, a*} of O's and l's such that 

Sq 2 a* =sq2[ ^,a*]. hclass [£, q + 2] . 

Then, we construct the matrix corresponding to Sq 2 H q K with the command 

matrixSq2 [£, q] 

=Table [sq2 U, hclass L£, q~\ [[«]]] 1, Length [hclass [£, q~\ ] } 

After this, we compute 

NullSpace [matrixSq2 l£, q~\ , Modulus— > 2] . hclass [£, q~\ 

in order to obtain a base of N q K . 

An example of the computation of Adem secondary cohomology operation 
using our algorithm is the following. Let if be a simplicial complex whose set 
of maximal simplices is 

{ (1, 3, 7), (3, 4, 7), (1, 4, 7), (1, 2, 8), (2, 3, 8), (1, 3, 8), 

(4, 5, 9), (4, 6, 9), (5, 6, 9), (3, 4, 10), (3, 6, 10), (4, 6, 10), 
(1, 2, 3, 4, 5, 6), (1, 2, 3, 4, 5, 11), (1, 2, 3, 4, 6, 11), 
(1,2,3, 5, 6, 11), (1,2,4, 5, 6, 11), (1,3, 4, 5,6, 11), (2,3, 4, 5,6, 11) } 
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We first compute the chain contraction to the homology: 

{ {(1),(1),0},{(2),(1),(1,2)},{(3),(1),(1,3)}, 

{(4), (1), (1,3) + (3, 4)}, {(5), (1), (1,3) + (3,4) + (4,5)}, 
{(6), (1), (1, 3) + (3, 4) + (4, 6)}, {(7), (1), (1, 7)}, {(8), (1), (1, 8)}, 
{(9), (1), (1, 3) + (3, 4) + (4, 9)}, {(10), (1), (1, 3) + (3, 10)}, 
{(11), (1), (1,11)}, {(1,4), 0,(1,3, 7) + (1,4, 7) + (3,4, 7)}, 
{(1,5), 0, (1,3,7) + (1,4,5) + (1,4,7) + (3,4,7)}, 
{(1, 6), 0, (1, 3, 6) + (3, 4, 10) + (3, 6, 10) + (4, 6, 10)}, 
{(2,3),0,(l,2,8) + (l,3,8) + (2,3,8)}, 
{(2, 4), 0, (1, 2, 4) + (1, 3, 7) + (1, 4, 7) + (3, 4, 7)}, 
{(2,5),0, (1,2,5) + (1,3,7) + (1,4,5) + (1,4,7) + (3,4,7)}, 
{(2, 6), 0, (1, 2, 6) + (1, 3, 6) + (3, 4, 10) + (3, 6, 10) + (4, 6, 10)}, 
{(2,8),0,(1,2,8)},{(2, 11), 0,(1,2,11)}, 
{(3,5),0, (1,3,5) + (1,3,7) + (1,4,5) + (1,4,7) + (3,4,7)}, 
{(3, 6), 0, (3, 4, 10) + (3, 6, 10) + (4, 6, 10)}, {(3, 7), 0, (1, 3, 7)}, 
{(3,8), 0,(1, 3,8)}, {(3, 11), 0,(1, 3,11)}, 
{(4, 7), 0, (1, 3, 7) + (3, 4, 7)}, {(4, 10), 0, (3, 4, 10)}, 
{(4, 11), 0, (1,3,7) + (1,4,7) + (1,4,11) + (3,4,7)}, 
{(5, 6), 0, (4, 5, 9) + (4, 6, 9) + (5, 6, 9)}, {(5, 9), 0, (4, 5, 9)}, 
{(5,11),0, (1,3,7) + (1,4,5) + (1,4,7) + (1,5,11) + (3,4,7)}, 
{(6, 9), 0, (4, 6, 9)}, {(6, 10), 0, (3, 4, 10) + (4, 6, 10)}, 
{(6, 11), 0, (1, 3, 6) + (1, 6, 11) + (3, 4, 10) + (3, 6, 10) + (4, 6, 10)}, 
{(1,2,3), (1,2,3), 0},{(1,3,4), (1,3, 4), 0}, 
{(1,4,6), (1,4, 6), 0},{(1,5,6), (1,5,6), 0}, 
{(2,3, 4), (1,2, 3) + (1,3,4), (1,2, 3,4)}, {(2, 3, 5), (1,2, 3), (1,2, 3, 5)}, 
{(2,3, 6), (1,2,3), (1,2, 3,6)}, {(2, 3, 11), (1,2,3), (1,2,3, 11)}, 
{(2,4, 5), 0,(1,2,4, 5)}, {(2,4, 6), (1,4, 6), (1,2, 4, 6)}, 
{(2,4, 11), 0,(1,2,4, 11)}, {(2, 5, 6), (1,5, 6), (1,2, 5, 6)}, 
{(2, 5, 11), 0,(1,2, 5, 11)}, {(2, 6, 11), 0,(1,2, 6, 11)}, 
{(3,4, 5), (1,3,4), (1,3, 4, 5)}, {(3,4, 6), (1,3,4) +(1,4, 6), (1,3,4, 6)}, 
{(3,4, 11), (1,3,4), (1,3, 4, 11)}, {(3, 5,6), (1,5,6), (1,3, 5,6)}, 
{(3, 5, 11), 0,(1,3, 5, 11)}, {(3,6, 11), 0,(1,3, 6, 11)}, 
{(4, 5,6), (1,4, 6) + (1,5, 6), (1,4, 5, 6)}, {(4, 5, 11), 0,(1,4, 5, 11)}, 
{(4, 6, 11), (1,4, 6), (1,4, 6, 11)}, {(5, 6, 11), (1,5, 6), (1,5, 6, 11)}, 
{(2,3, 4, 5), 0,(1,2,3, 4, 5)}, {(2,3, 4, 6), 0,(1,2, 3, 4, 6)}, 
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{(2, 3,4, 11), 0,(1,2, 3,4, 11)}, {(2, 3, 5,6), 0,(1,2,3, 5, 6)}, 
{(2, 3, 5, 11), 0,(1,2, 3, 5, 11)}, {(2,3, 6, 11), 0,(1,2,3, 6, 11)}, 

{(2, 4, 5, 6), 0, (1, 2, 4, 5, 6)}, {(2, 4, 5, 11), 0, (1, 2, 4, 5, 11)}, 
{(2, 4, 6, 11), 0,(1,2,4, 6, 11)}, {(2, 5, 6, 11), 0,(1,2, 5, 6, 11)}, 
{(3, 4, 5, 6), 0,(1,3, 4, 5, 6)}, {(3, 4, 5, 11), 0,(1,3, 4, 5, 11)}, 
{(3, 4, 6, 11), 0,(1,3, 4, 6, 11)}, {(3, 5,6, 11), 0,(1,3, 5, 6, 11)}, 
{(4, 5, 6, 11), 0,(1,4, 5, 6, 11)}, {(2, 3, 4, 5, 6), 0,(1,2,3, 4, 5,6)}, 
{(2, 3,4, 5, 11), 0,(1,2,3, 4, 5, 11)}, {(2,3, 4, 6, 11), 0,(1,2, 3, 4, 6, 11)}, 
{(2, 3, 5, 6, 11), 0,(1,2,3, 5,6, 11)}, {(2,4, 5, 6, 11), 0,(1,2,4, 5,6, 11)}, 
{(3, 4, 5, 6, 11), 0,(1,3,4, 5, 6, 11)}, {(2,3, 4, 5,6, 11), (2,3, 4, 5, 6, 11), 0} }. 

Notice that if a simplex of K doesn't appear in this list, it is because its image 
under / and cf) is null. The representative cycle of every homology class is: 

.9(1) = (1) 

ff (l,2,3) = (1,2,3) + (1,2,8) + (1,3,8] + (2,3,8) 
ff (l,3,4) = (1,3,4) + (1,3,7) + (1,4,7) + (3,4,7) 

ff (l, 4, 6) = (1,3,4) + (1, 3, 6) + (1, 4, 6) + (3, 4, 10) + (3, 6, 10) + (4, 6, 10) 
3(1, 5, 6) = (1, 4, 5) + (1, 4, 6) + (1, 5, 6) + (4, 5, 9) + (4, 6, 9) + (5, 6, 9) 
3(2, 3, 4, 5, 6, 11) = (1, 2, 3, 4, 5, 6) + (1, 2, 3, 4, 5, 11) + (1, 2, 3, 4, 6, 11) 
+ (1,2, 3, 5, 6, 11) + (1,2,4, 5,6, 11) + (1,3, 4, 5, 6, 11) + (2, 3,4, 5,6, 11). 

A base of the kernel of Sq 2 H 2 K is: 

{(1,2,3)*, (1,3,4)*, (1,4,6)*, (1,5,6)*}. 

Now, given an element a of this kernel, we first have to compute c = g*a. Let 
us study a concrete example with all the details. Let us take a = (1,2,3)* + 
(1,5,6)*. Then 

c = g*a = af = (1, 2, 3)* + (1, 5, 6)* + (2, 3, 4)* + (2, 3, 5)* + (2, 3, 6)* 

+ (2, 3, 11)* + (2, 5, 6)* + (3, 5, 6)* + (4, 5, 6)* + (5, 6, 11)* . 

We now compute the cochains of the 3rd step of the algorithm for computing 

6b = c^c = (1,2,3,5,6)* + (2,3,4,5,6)* + (2,3,5,6,11)* 

b = (c — c)<j> = (2,3,5,6)* 

Then, we have that b — i 6 = and b - 2 Sb = 0. On the other hand, Si] = 
c ^— "i c = therefore i] 8rj = 0. We thus get, 

w = f*(E 3 c 4 ) = (Ezc A )g = (1, 2, 3, 4, 5, 6)* 5 = (2, 3, 4, 5, 6, 11)* . 
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Therefore, tf 2 «l,2,3)* + (1,5,6)*) = (2,3,4,5,6,11)*. Finally, observe that 
since there are no classes of cohomology of dimension 3, then (2, 3, 4, 5, 6, 11)* 
ImSq 2 H 3 K. 

6 Some Comments 

All these results can be given in a more general framework working not nec- 
essarily with finite simplicial complexes. Nevertheless, a contraction from the 
(co)chain complex associated to the simplicial complex to its (co)homology must 
exist in order to develop the method. 

In this paper, the ground ring is Z 2 for simplicity, but the same process 
can be done working with any held as the ground ring. For example, let Z p (p 
being a prime) be the group of coefficients. From the combinatorial formulae 
for the reduced pth powers [Ste47, SE62] at cochain level in terms of face 
operators established in [GR99, GonOO] and the algorithm for computing the 
chain contraction (/, g, (f>) from C*(K; Z p ) to H^,(K; Z p ), Steenrod cohomology 
operations can effectively be computed. Let a* G H q (K;Z p ), for calculating 
the cohomology class Vi{a*) with a* E H 9 (K;Z p ), we only have to compute 
Pi{af)g. 

Finally, in order to obtain the image of any cohomology operation at cochain 
level over a representative cocycle using our formulae, we have to compute them 
on a base of C* (K) in the desired dimension. A way of decreasing the complexity 
of this is to do a "topological" thinning of the simplicial complex K in order 
to obtain a thinned simplicial subcomplex M top K of K (such that there exists 
a chain contraction from C*K to C*(M iop K)). Two examples of thinning in 
this way are edge contractions (example (a)) and simplicial collapses (example 
(b)). Therefore, we can apply our machinery to compute cohomology operations 
in the thinned simplicial complex M top K and then, the results can be easily 
interpreted in the "big" simplicial complex K via composition of contractions. 
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